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Although animals eliminate apoptotic cells using 
macrophages, plants use cell corpses throughout 
development and disassemble cells in a cell- 
autonomous manner by vacuolar cell death. During vacu- 
olar cell death, lytic vacuoles gradually engulf and digest 
the cytoplasmic content. On the other hand, acute stress 
triggers an alternative cell death, necrosis, which is char- 
acterized by mitochondrial dysfunction, early rupture of 
the plasma membrane, and disordered cell disassembly. 
How both types of cell death are regulated remains ob- 
scure. In this paper, we show that vacuolar death in the 



embryo suspensor of Norway spruce requires autophagy. 
In turn, activation of autophagy lies downstream of meta- 
caspase mcll-Pa, a key protease essential for suspensor 
cell death. Genetic suppression of the metacaspase- 
autophagy pathway induced a switch from vacuolar to 
necrotic death, resulting in failure of suspensor differentiation 
and embryonic arrest. Our results establish metacaspase- 
dependent autophagy as a bona fide mechanism that is 
responsible for cell disassembly during vacuolar cell death 
and for inhibition of necrosis. 



Introduction 

Programmed cell death (PCD) is indispensable for animal and 
plant development, but the mechanisms of PCD differ between 
the two kingdoms. Plants lack apoptosis that involves cell frag- 
mentation into discrete bodies and their heterophagic removal, 
owing to the presence of cell walls and lack of phagocytosis 
(Beers, 1997; Jones, 2001; Lam, 2004). Furthermore, plant ge- 
nomes lack the core apoptotic regulators, such as Bcl-2 family 
proteins and caspases (Koonin and Aravind, 2002). 

Although molecular regulation of plant PCD remains 
poorly understood, most cases of plant cell death can be divided 
into two classes with distinct kinetics and morphology: vacuolar 
cell death and necrosis (van Doom et al., 2011). Vacuolar cell 
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death is a slow process whereby growing lytic vacuoles gradu- 
ally digest entire or most of the contents of terminally differen- 
tiated cells excluding cell walls. This cell death is indispensable 
for plant development, playing an instrumental role in the for- 
mation of conduits of water, nutrients, and hormones (the em- 
bryo suspensor and the vascular system) and secretory structures 
(e.g., laticifers; Beers and McDowell, 2001; Bozhkov et al., 
2005a; van Doom and Woltering, 2005; Bollhoner et al., 2012). 
We have shown that execution of vacuolar cell death in Norway 
spruce (Picea abies) embryo suspensor requires type II meta- 
caspase mcll-Pa, a cysteine protease distantly related to cas- 
pases (Suarez et al., 2004; Bozhkov et al., 2005b; Tsiatsiani 
et al., 201 1). However, mcll-Pa is a nucleocytoplasmic protein, 
and its proteolytic action alone is not sufficient to account for 
progressive growth of lytic vacuoles in the course of execution 
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of PCD. Identification of such a mechanism would significantly 
advance our understanding of the vacuolar cell death process. 

In contrast to vacuolar cell death, necrosis is a rapid mode 
of plant cell death and, like in animal systems, involves mito- 
chondrial dysfunction and early loss of plasma membrane in- 
tegrity (van Doom et al., 201 1; Galluzzi et al., 2012). However, 
although animal necrotic cells swell and burst, their plant coun- 
terparts retain undigested cytoplasmic materials surrounded by 
cell walls. Necrosis is a typical response of plants to abiotic 
stress and pathogen attack (Heath, 2000; Coll et al., 2011; 
van Doom et al., 2011). Thus, vacuolar and necrotic modes of 
cell death not only differ morphologically but are also impli- 
cated in distinct aspects of plant biology: development and 
stress response, respectively. An intriguing question is what de- 
termines the mode of cell death in plants. 

Here, we show that autophagy is responsible for cell self- 
disassembly during vacuolar cell death in P. abies embryos. 
Activation of autophagy requires metacaspase mcll-Pa and 
deficiency of either component switches the mode of cell death 
from vacuolar to necrotic. These findings provide a mechanistic 
explanation for morphological differences between two major 
classes of cell death in plants. 

Results and discussion 

Vacuolar cell death in the embryo 
suspensor is associated with 
enhanced autophagy 

In somatic embryogenesis of P. abies, the formation of early 
embryos is induced by withdrawal of growth factors (GFs), 
auxin and cytokinin, whereas further development to the cotyle- 
donary stage requires exogenous abscisic acid (ABA; Fig. 1 A; 
Bozhkov et al., 2005a). An early P. abies embryo is composed 
of a proliferating embryonal mass (EM) that will eventually 
form a cotyledonary embryo and terminally differentiated sus- 
pensor, which is gradually eliminated before the cotyledonary 
stage. Although Arabidopsis thaliana embryos have minute sus- 
pensors of seven cells, the suspensors in P. abies and most other 
gymnosperms are several millimeters long and composed of 
many cells (Fig. 1 B; Singh, 1978). In addition, suspensors of 
P. abies consist of several tiers of elongated cells at successive 
stages of cell disassembly, providing an excellent paradigm 
for studying vacuolar PCD (Bozhkov et al., 2005a; van Doom 
et al., 2011). 

We obtained three lines of evidence that vacuolar PCD in 
the P. abies suspensor is associated with increased autophagic 
activity. First, transmission electron microscopy (TEM) re- 
vealed accumulation of autophagic bodies in the vacuoles of 
suspensor cells upon inhibition of vacuolar acidification using 
concanamycin A (ConA; Fig. 2 A) as well as increased amounts 
of double membrane-bound autophagosomes in the cytoplasm 
of suspensor cells as compared with EM cells (Fig. 2, A-C; 
Filonova et al., 2000). Second, transgenic mRFP-Atg8 lines 
showed cytoplasmic localization of mRFP-Atg8 in the EM 
cells and punctate localization in the suspensor cells (Fig. 2 D; 
Klionsky et al., 2012). Simultaneous measurement of fluores- 
cein diacetate (FDA) staining intensity, cell length, and amount 




Figure 1 . Embryo development in P. abies. (A) A schematic model of 
somatic embryogenesis. Cell proliferation is stimulated by GF, auxin, and 
cytokinin. Early embryogenesis is induced by withdrawal of GF. An early 
embryo is composed of embryonal mass (EM) and a suspensor. Devel- 
opment to the cotyledonary stage is promoted by abscisic acid (ABA). 
Suspensors are eliminated before the cotyledonary stage. Dying or dead 
cells are colored in blue. (B) Embryos of P. abies and A. thaliana (inset; 
dashed lines indicate contour of suspensor) at the developmental stage 
before formation of cotyledons stained with fluorescein diacetate (FDA; 
green), DAPI (blue), and FM4-64 (red). The lack of FDA staining in the 
suspensor denotes the loss of cell viability. Note the giant size, as well as 
the higher suspensor-to-EM size ratio, for P. abies embryo as compared 
with the A. thaliana embryo. Bars, 50 urn. 

of mRFP-Atg8 puncta per cell area in the EM and suspensor 
cells confirmed that progression of vacuolar PCD in the suspen- 
sors correlates with cell elongation and enhanced autophagy 
(Fig. 2 E). Finally, abrogation of autophagic flux by ConA led 
to dramatic increase in the levels of autophagic target proteins 
Atg8 and NBR1 (Fig. 2 F; Svenning et al., 2011; Klionsky et al., 
2012;Minina et al, 2013b). 

Silencing of AT (35 or AT OS abrogates 
embryogenesis and phenocopies 
mcll-Pa depletion 

The vast body of genetic data demonstrates conservation of 
the major functional groups of Atg (autophagy related) pro- 
teins in plants (Doelling et al., 2002; Hanaoka et al., 2002; 
Liu and Bassham, 2012). We addressed the requirement of au- 
tophagy for vacuolar PCD in the suspensor by silencing two 
single-copy ATG genes, ATG5 and ATG6, by RNAi (Figs. SI 
and S2). Silencing either gene abrogated suspensor formation 
and apical-basal patterning (Fig. 3 A), the phenotype being 
similar to lines deficient for mcll-Pa, the only type II metacas- 
pase expressed in the P. abies embryos (Figs. 3 A, SI, and S2; 
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Figure 2. Enhanced autophagy in the P. abies embryo suspensor. (A) Assessment of autophagic flux in the EM and suspensor cells using ConA treat- 
ment. Arrows denote autophagic bodies. Insets depict autophagosome docking to the vacuole (—ConA) and autophagic body (+ConA). N, nucleus; 
V, vacuole; asterisks, cell wall. Bars, 2 urn. (B) Typical autophagosome in the suspensor cell. Bar, 0.2 urn. (C) Number of autophagosomes in the cytoplasm 
estimated from the micrographs of EM and suspensor (S) cells. Data represent the means ± SEM from three independent experiments, each including at 
least four cells per cell type. (D) mRFP-Atg8 accumulates in puncta in the suspensor cells but remains cytoplasmic in the EM cells upon ConA treatment. 
Bars, 50 pm. (E) Correlation between anisotropic cell expansion (cell length), progression of cell death (mean FDA intensity; green symbols and green 
power trendline, R 2 = 0.5402), and accumulation of autophagosomes (number of mRFP-Atg8-positive puncta; red symbols and red polynomial trendline, 
R 2 = 0.8202) in the embryos. RFU, relative fluorescence units. (F) Accumulation of NBR1 , mRFP-Atg8, and Atg8 upon ConA treatment detected by West- 
ern blot analysis of total protein extracts from the embryos. Actin was used as a reference control. Graph represents relative protein amounts shown as 
means ± SEM from three independent measurements. For each protein, the integrated band intensities were first normalized to corresponding intensities 
of Coomassie staining (loading) and then to the sample without ConA. The amounts of all three proteins were significantly increased upon ConA treat- 
ment (P < 0.001; Student's ftest). 
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Figure 3. Autophagy and metacaspase are essential for embryo patterning. (A-D) Knockdown of ATG5, ATG6, or mcll-Pa results in suspensor-deficient 
phenotype (A) and inhibits anisotropic growth of dying cells (B), early (C), and late (D) embryogenesis. (A) The cultures were grown without GF for 7 d and 
stained with Evans blue to detect dead cells. Bars, 20 urn. (B) Length of Evans blue-positive cells measured in the same samples. Data represent the means ± 
SEM from three independent experiments, each including at least 200 Evans blue-positive cells from 1 0 embryos per line. (C) Frequency of aberrant early 
embryos lacking suspensors in the same samples as in A. Data represent the means ± SEM from three independent experiments, each including >100 
embryos per line. (D] The yield of cotyledonary embryos after ABA treatment. Data represent the means ± SEM from four independent experiments, each 
including eight samples per line. FW, fresh weight. The mean values for RNAi lines in B-D significantly differ from control (P < 0.001 ; Dunnett's test). 



Suarez et al., 2004). The failure of normal suspensor formation 
in Atg5-, Atg6- and mcII-Pa-deficient lines was accompanied 
by reduced anisotropic expansion of dying cells as estimated 
by measuring the length of Evans blue-positive cells (Fig. 3 B). 
Genetic inhibition of autophagy or metacaspase increased the 
frequency of aberrant early embryos lacking suspensors in sev- 
eral independent lines by four- to sixfold (Fig. 3 C). Accord- 
ingly, the efficiency of subsequent embryo development to 
the cotyledonary stage was dramatically reduced (Fig. 3 D; 
Suarez et al., 2004). 



Deficiency of either autophagy or 
metacaspase switches the mode of 
cell death from vacuolar to necrotic 

Because autophagy or metacaspase deficiency did not prevent 
cell death upon withdrawal of GF (Fig. 3 A), we compared the 
morphology of cell death in RNAi and control lines using TEM. 
In control lines, dying suspensor cells maintained intact plasma 
membranes and turgid protoplasts until almost all of the cyto- 
plasm was digested by lytic vacuoles, whereas dying cells in 
RNAi lines contained shrunken and largely undigested protoplasts 
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with ruptured plasma membranes (Fig. 4, A-D). This morphol- 
ogy was reminiscent of necrosis (Majno and Joris, 1995; van 
Doom et al., 2011). The reduced number of morphologically 
intact mitochondria in the cytoplasm of dying cells from the 
RNAi lines as compared with both meristematic and suspensor 
cells (Fig. 4 E), correlated with the accumulation of swollen 
mitochondria (Fig. 4, F and G). Consistent with the altered mor- 
phology of mitochondria, the RNAi lines displayed a drop in 
intracellular ATP content (Fig. 4 H) and respiratory decline 
(Fig. 4 I; Eguchi et al., 1997; Bal-Price and Brown, 2000). 
These metabolic changes indicate mitochondrial dysfunction 
and, in combination with the loss of plasma membrane integ- 
rity, demonstrate that the cell death observed in autophagy- or 
metacaspase-deficient lines is necrosis. 

The main reason for the incomplete clearance of cellular 
contents in necrotic cells (Fig. 4 C) was compromised growth 
of lytic vacuoles, as evidenced from the measurements of the 
ratios of vacuolar to the whole protoplast areas on the electron 
micrographs (Fig. 4 J). The failure of progressive vacuolization 
of necrotic cells correlated with decreased autophagic activity 
in the cytoplasm (Fig. 4 K) compared with the suspensor cells 
from control lines in which enhanced autophagy sustained vac- 
uole-mediated cell clearance (Fig. 4, K and L). 

To quantify necrotic death in different lines, embryogenic 
cultures were stained with FM4-64 to distinguish between in- 
tact and shrunken protoplasts and FDA to assess cell viability. 
Combination of both dyes provides robust diagnostic tool for 
identification of necrotic cells (Fig. 4 M; Minina et al., 2013a). 
Although the frequency of necrotic cells (as distinguished by 
collapsed protoplast and no FDA staining) was negligible in the 
control lines, it increased to 60% in the individual RNAi lines 
(Fig. 4 N). No statistical differences were found between lines 
silenced for ATG5, ATG6, or mcll-Pa. Collectively, our data 
demonstrate that autophagy and metacaspase mcll-Pa are re- 
quired for both progression of vacuolar cell death and suppres- 
sion of necrosis. 

Metacaspase lies upstream of autophagy 
in the vacuolar cell death pathway 

We have previously shown that the physiological effect of 
mcll-Pa in the control of suspensor PCD requires its proteo- 
lytic activity (Bozhkov et al., 2005b). To address whether 
metacaspase lies up- or downstream of autophagy in the vacu- 
olar cell death pathway, we first compared proteolytic activity 
of mcll-Pa in the protein extracts from control, ATG5 RNAi, 
and ATG6 RNAi lines using three independent approaches: 
cleavage of peptidic substrate FESR-7-amino-4-methylcoumarin 
(AMC; Fig. 5 A), mcll-Pa autoprocessing, and cleavage of Tudor 
staphylococcal nuclease (TSN), a target of mcll-Pa (Fig. S3; 
Bozhkov et al., 2005b; Sundstrom et al., 2009). Because there 
was no difference in the levels of metacaspase activity observed 
between autophagy-deficient and control lines using all three 
approaches, we rejected the possibility that autophagy lies up- 
stream of mcll-Pa in the pathway. 

To test an alternative scenario in which mcll-Pa is upstream 
of autophagy, we expressed mRFP-Atg8 in the wild-type (con- 
trol) and RNAi lines and quantified the amount of autophagosomes 



per volume of dying cells by confocal microscopy after ConA 
treatment (Fig. 5, B and C). As shown in Fig. 5 C, silencing 
of mcll-Pa led to a dramatic decrease in the number of autopha- 
gosomes in the dying cells, to levels found in Atg5- and Atg6- 
deficient lines. Suppression of autophagy under mcll-Pa deficiency 
was further confirmed by a significant decrease in the efficiency 
of the degradation of NBR1 (Fig. 5 D). Together, these data 
establish metacaspase as an upstream mediator of autophagy in 
the vacuolar cell death pathway (Fig. 5 E). 

Conclusions 

Our combined data demonstrate a cell-autonomous require- 
ment of autophagy and mcll-Pa for the execution of vacuolar 
cell death in the P. abies embryo suspensor, which in turn is 
essential for embryogenesis. Critical dependence of vacuolar 
death on autophagy and type II metacaspases seems to be in- 
dependent of the experimental system because progression of 
xylem PCD (classic example of vacuolar cell death; van Doom 
et al., 2011) in A. thaliana was shown to require autophagy 
(Kwon et al., 2010) and type II metacaspase AtMC9 (Bollhoner 
et al., 2013), although these two components have not been 
linked together in a pathway. The contribution of vacuolar cell 
death to developmental processes in different systems can 
vary. For example, the knockout mutants of ATG genes in 
A. thaliana are fertile and do not reveal massive embryo abortion 
caused by autophagy deficiency in the suspensor. The differ- 
ence between A. thaliana and P. abies can be explained by the 
minute contribution of the seven-cell suspensor to the A. thali- 
ana embryo in comparison to the giant suspensors of P. abies 
that exceed the EMs (Fig. 1 B). Furthermore, although the sus- 
pensor cells in P. abies do not divide and are subjected to PCD 
as soon as they are formed (Smertenko et al., 2003; Bozhkov 
et al., 2005a), the suspensor cells in A. thaliana undergo sev- 
eral rounds of divisions before PCD is initiated during transi- 
tion from the late globular to the heart stages (Mansfield and 
Briarty, 1991; Kawashima and Goldberg, 2010; Wendrich and 
Weijers, 2013). This implies that in addition to size differ- 
ences, the suspensor in P. abies embryos is composed entirely 
of dying and dead cells regardless of the developmental stage, 
whereas the fate of suspensor cells in A. thaliana is determined 
by the stage of embryo development. 

Our work strengthens the idea about the evolutionarily 
conserved and context-independent role of autophagy in the 
suppression of necrosis and establishes a new paradigm for 
unraveling molecular mechanisms underlying relationships be- 
tween two major classes of cell death in plants. In vacuolar cell 
death, autophagy is both an executioner of slow yet complete 
cell clearance and a protector against rapid cell death, necrosis, 
which leaves behind a largely unprocessed cell corpse. Autophagy 
was also shown to counteract the spread of necrotic lesions in 
plant leaves infected with microbial or fungal pathogens (Liu 
et al., 2005; Hofius et al., 2009; Yoshimoto et al., 2009; Lenz 
et al., 201 1). Finally, in mammalian systems and Dictyostelium 
discoideum, autophagy differentially controls distinct modes of 
cell death and usually suppresses necrosis (Kosta et al., 2004; 
Degenhardt et al., 2006; Chen et al., 2008; Kroemer and Levine, 
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Figure 4. Deficiency of either autophagy or metacaspase switches the mode of cell death from vacuolar to necrotic. (A-C) Cell morphology in control and 
RNAi lines (exemplified by an ATG6 RNAi line). N, nucleus; V, vacuole; asterisks, cell wall. (A) A meristematic cell from the EM of the control line. 
(B) Vacuolated suspensor cell with intact plasma membrane and turgid protoplast. (C) Dead cells from an RNAi line with shrunken protoplasts (double 
arrows denote detachment of plasma membranes from the cell walls) and ruptured plasma membranes (inset, arrow). (D) Frequency of cells with com- 
promised plasma membranes assessed by SYTOX orange staining of meristematic (M; control line), suspensor (S; control line), and necrotic (N; RNAi 
line) cells. Data represent the means ± SEM from assessing all cells in >10 embryogenic structures per line. (E-G) Mitochondria in control and RNAi lines 
(exemplified by an ATG6 RNAi line). (E) Number of intact mitochondria per micrometer squared of the cytoplasm counted on the micrographs of meriste- 
matic, suspensor, and necrotic cells. Data represent the means ± SEM from three independent experiments, each including at least four cells per cell type. 
(F) Intact mitochondria (arrows) in a suspensor cell from the control line. (G) Swollen and degraded mitochondria (arrows) in a cell from an RNAi line. 
(H and I) Intracellular ATP content (H) and oxygen consumption (I). Data represent the means ± SEM from three independent experiments. *, P < 0.05; 
**, P < 0.01 (vs. control, Dunnett's test). FW, fresh weight. (J and K) Cell vacuolization [I] and number of autophagosomes in the cytoplasm (K) estimated from 
the micrographs of meristematic (control line), suspensor (control line), and necrotic [ATG6 RNAi line) cells. Data represent the means ± SEM from 
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Figure 5. Metacaspase lies upstream of autophagy in the vacuolar cell death pathway. (A) Proteolytic activity of mcll-Pa in control and RNAi lines mea- 
sured using the FESR-AMC peptide. Data represent the means ± SEM from three independent experiments. The mean values were normalized to the control 
line. *, P < 0.05; **, P < 0.01 (vs. control, Dunnett's test]. (B) Localization of mRFP-Atg8 in dying cells from control and RNAi lines with and without ConA 
treatment. Bars, 50 urn. DIC, differential interference contrast microscopy. (C) Quantification of autophagosomes. Data represent the means ± SEM number 
of mRFP-Atg8-positive puncta counted for >10 cell volumes, 10 urn 3 each, per line. **, P < 0.01 ; ***, P < 0.001 (vs. control, Dunnett's test]. (D) Western 
blot analysis of NBR1 degradation. Numbers shown below the blots correspond to the relative amount of NBR1 . The integrated band intensities were first 
normalized to corresponding intensities of Coomassie staining (loading) and then to the control. (E) Schematic pathway showing the relationship between 
autophagy and metacaspase mcll-Pa in the regulation of vacuolar and necrotic modes of cell death. 



2008; Yuan and Kroemer, 2010; Shen et al., 2012). Interestingly, 
although autophagy-dependent vacuolar cell death in plants 
determines the cell fate (e.g., differentiation of the suspensor 
cells) and cannot be substituted by necrosis, the switch from 
autophagy-dependent to -independent cell death in mammals re- 
sults in no major developmental defects (Golstein and Kroemer, 
2005). This indicates that the role of autophagic cell death in 
development is not conserved throughout the evolution and can 
be redundant. 



Our study shows that autophagy is an integral part of the 
suspensor cell death in P. abies; however, instead of being the 
effector mechanism of cell death, it performs cytoprotective 
role and orchestrates slow kinetics and morphological features 
of the vacuolar cell death. This demonstrates that the suspen- 
sor cells die with autophagy but not by autophagy (Kroemer 
and Levine, 2008). The primary and so far unknown cell death 
induction signal, as well as point of no return in the cell death 
pathway, lies apically to both autophagy and its upstream 



three independent experiments, each including at least four cells per cell type. (L] Autophagosomes (open arrows) in the cytoplasm and autophagic bodies 
(closed arrows) in the vacuole of the embryo suspensor cell from a control line. (M and N) Quantification of necrotic death in control and RNAi lines after 
staining with FDA and FM4-64. (M) Staining patterns of EM and suspensor cells (control line) and necrotic cells (RNAi lines). DIC, differential interference 
contrast microscopy. (N) Frequency of necrotic cells. Data represent the means ± SEM from two independent experiments, each including two lines per 
genotype and >50 cells per line (P < 0.001 ; vs. control, Dunnett's test). Different letters in D, E, J, and K indicate statistically different means (P < 0.001 ; 
Student's ftest). Bars: (A-C, main image) 5 urn; (C, inset) 0.5 pm; (F, G, and L) 2 urn; (M) 10 urn. 
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regulator metacaspase mcll-Pa (Fig. 5 E). A similar hierarchy 
was observed in PCD during Drosophila meJanogaster oogene- 
sis, in which caspases act upstream to autophagy (Hou et al., 
2008). Further understanding of proteolytic processes mediated 
by evolutionarily related caspases and metacaspases is required 
to identify the molecular mechanisms activating autophagy dur- 
ing execution of PCD in animals and plants. 

Materials and methods 

Phylogenetic analysis of Atg5, Atg6, and metacaspase homologues 

A search for proteins containing PF00656 (CI 4 peptidase), PF04106 
(Apg5 domain), or PF041 1 1 (Apg6 domain) was performed in Pfam data- 
base for species enlisted in the Table SI . An alternative search for proteins 
containing cl00042 (CASc superfamily) or cl04402 (APG5 superfamily) 
was performed using CDART (Conserved Domain Architecture Retrieval 
Tool) tool (Geer et al., 2002). Obtained hits were used for multiple align- 
ments in ClustalW to confirm the presence of p20 and pi 0 domains in meta- 
caspase homologues, Beclin-1 domain in Atg6 homologues, and Apg5 
domain in Atg5 homologues. Estimated conserved sequence elements from 
the alignments were used as a query for additional search of homologues 
from gymnosperm species using TBtASTN tool in NCBI and ConGenlE data- 
bases. To confirm the absence of additional homologues, verified sequences 
for each individual species were cross-checked in proteome, transcriptome, 
and genome databases available at Phytozome, Ensembl, and ConGenlE. 
Multiple alignments of verified sequences were performed in ClustalW. Un- 
rooted trees were constructed using the neighbor-joining method (Saitou and 
Nei, 1987) using the yeast homologue as an out group. Topology of each 
tree was confirmed by bootstrap analysis (2,000 repeats). 

Cloning of hairpin and overexpression constructs 

Full-length cDNA sequences of P. abies ATG5 (available from GenBank 
under accession no. HE793992.1) and ATG6 (GenBank accession no. 
HG326275) were obtained by rapid amplification of cDNA ends PCR 
using specific primers designed based on comparison of full-length or 
partial sequences of homologues from Oryza sativa, A. thaliana, Picea 
glauca, and Picea sitchensis. To produce the hairpin of ATG5, the correspond- 
ing part of ATG5 cDNA was amplified using primers hpATG5-P-sense/ 
ATG5-reverse and hpATG5-M-antisense/ATG5-reverse (the sequences of 
all primers are listed in Table S2) and cut with Bglll-Ascl and Pacl-BamHI 
(Thermo Fisher Scientific), respectively, producing 1 38-635-nucleoside tri- 
phosphate (ntp) and 138-847-ntp fragments of ATG5 coding DNA se- 
quence (CDS). The fragments were ligated with T4 DNA ligase (Thermo 
Fisher Scientific) into a modified pAHC25 vector cut with Pacl-Ascl. The 
full-length hairpin was then amplified with primers attBl -pAHC25-forward/ 
attB2-pAHC25-reverse and recombined into pDONR/Zeomycin (Zeo) vec- 
tor (Invitrogen) followed by tR recombination into the pMDC32 vector 
(Curtis and Grossniklaus, 2003). To produce the hairpin of ATG6, the 
corresponding part of ATG6 cDNA was amplified using primers Beclin- 
Ascl-sense/attB2-hpATG6 and attBl -hpATG6/Beclin-Ascl-Ant. Both PCR 
products were then cut with the AscI enzyme (Thermo Fisher Scientific) pro- 
ducing 689-1, 298-ntp and 826-1 ,298-ntp fragments of ATG6 CDS, 
which were then ligated into pAHC25 and cloned into the pDONR/Zeo 
vector followed by recombination into pMDC32 vector as described for 
ATG5. To produce the hairpin of mcll-Pa, fragments of mcll-Pa cDNA (GenBank 
accession no. AJ534970.2) were amplified using primers attBl -hpmcll-Pa/ 
as-hpmcll-Pa-EcoRI and s-hpmcll-Pa-EcoRI/attB2-hpmcll-Pa and cut with 
EcoRI (Thermo Fisher Scientific) producing 3 1 2-1 ,253-ntp and 537-1, 
253-ntp fragments of mcll-Pa CDS, respectively. The PCR products were 
then ligated, and the resulting construct was recombined into pDONR/Zeo 
vector followed by !R recombination into pMDC32 vector. Full-length CDS 
of ATG8 (Iclcomp43277_c0_seql; Nystedt et al., 2013) was amplified 
using primers attBl -ATG8-Pa/attB2-ATG8-Pa, cloned into pDONR/Zeo 
vector and then into a modified pMDC32 vector containing mRFP. The 
GUS gene was amplified from the pAHC25 plasmid (Christensen and 
Quail, 1996) using primers attBl-GUS-forward/attB2-GUS-forward and 
cloned into the pDONR/Zeo vector and then into the pMDC32 vector. 

P. abies embryogenesis system 

P. abies embryogenic cell lines were subcultured every 7 d by transferring 
~3 ml of the cell pellet into 100 ml of fresh half-strength LP medium and 
incubated at 22°C in the darkness on a shaker at 120 rpm as described 
previously (Filonova et al., 2000; Bozhkov et al., 2002; Minina et al., 



201 3a). Early embryogenesis was induced by omitting GF auxin and cyto- 
kinin from the liquid medium (— GF medium) for 7 d. To assess the effi- 
ciency of early embryogenesis, 5 ml of 7-d-old — GF culture was stained 
with Evans blue, and the percentage of polarized and aberrant embryos 
was counted. To assess the efficiency of late embryogenesis, 7-d-old — GF 
suspension cultures were plated on filter papers placed on the top of solidi- 
fied maturation medium containing ABA (Filonova et al., 2000) and incu- 
bated in the darkness at 22°C. Fully developed cotyledonary embryos 
were scored in week 4 and 6 of ABA treatment. 

Transgenic lines 

Cell lines expressing mRFP-Atg8 in ATG5 RNAi, ATG6 RNAi, mcll-Pa RNAi, or 
control (expressing GUS) background were established by agrotransformation 
of P. abies wild-type embryogenic cell line 06.28.05. The constructs were 
transformed into the Agrobacterium tumefaciens c58cl strain carrying addi- 
tional virulence plasmid pTOK47. Transformed bacteria were then grown 
overnight in YEP (yeast extract peptone) medium supplemented with 1 00 mg/ 
liter rifampicin, 5 mg/liter tetracycline, 100 mg/liter carbenicillin, and 
50 mg/liter kanamycin. Cells were centrifuged at 4,500 g for 5 min and re- 
suspended in infiltration buffer (10 mM MgCI 2 , 10 mM MES, pH 5.5, and 
150 pM acetosyringone). The OD600 of the cultures was adjusted to 10 with 
the infiltration buffer. A. tumefaciens transformed with mRFP-Atg8 was added 
in a 1 : 1 ratio to the agrobacterial cultures transformed with GUS or RNAi con- 
structs and incubated at room temperature for 1 h. 1 ml of the A. tumefaciens 
mixture was added to 5 ml of a cell pellet of 7-d-old suspension culture grown 
in the presence of GF. The volume of the mixture was adjusted to 10 ml with 
GF-containing medium (+GF medium), and acetosyringone was added to a 
final concentration at 1 50 uhA. After 4 h, the cells were plated on a filter paper 
placed on the top of solidified +GF medium and incubated at room tempera- 
ture in the darkness for 48 h. Then, filters were transferred on +GF plates con- 
taining 400 mg/ml timentin and 250 mg/ml cefotaxime and incubated under 
the same conditions for 5 d. Subsequently, filter papers were transferred onto 
fresh +GF plates containing 20 ug/ml hygromycin, 400 ug/ml timentin, and 
250 ug/ml cefotaxime and subcultured onto the same medium every week. 
The transgenic colli were picked from the plates after 3 wk. The presence of 
hairpins was verified by PCR with genomic DNA. The genomic DNA was ex- 
tracted using a modified hexadecyltrimethylammonium bromide method as 
previously described (Ciavatta et al., 2002) and used for PCR with primers 
attBl /PaATG5-forward, attBl /Beclin-Ascl-antisense, and attBl /as-hpmcll- 
Pa-EcoRI for ATG5 RNAi, ATG6 RNAi, and mcll-Pa RNAi lines. Absence of 
A. tumefaciens in the cultures was confirmed by performing PCR with picA- 
forward/picA-reverse and VirG-forward/VirG-reverse primers. GUS expres- 
sion was confirmed by a histochemical assay as previously described by 
Jefferson et al. (1987). mRFP-Atg8 expression was confirmed by detecting 
fluorescence using a confocal microscope (LSM 780; Carl Zeiss), with excita- 
tion of 561 nm and emission of 563-686 nm. The wild-type line and lines ex- 
pressing both GUS and mRFP-Atg8 were used as controls. 

Culture staining and confocal microscopy 

All imaging was performed at room temperature. Samples were mounted 
in aqueous solutions. 5 ml of 7-d-old — GF suspension cultures were incu- 
bated with 0.0025% (wt/vol) Evans blue solution for 1 0 min on the bench, 
transferred onto a 50-um nylon mesh, and washed twice with — GF me- 
dium. For quantification of normal and aberrant cell structures, the cultures 
were embedded in 6% (wt/vol) low melting temperature agarose. Pheno- 
types were assessed using an inverted microscope (5x Achrostigmat and 
lOx Plan-Neofluar objectives with numerical apertures 0.12 and 0.3, re- 
spectively; Axiovert 10; Carl Zeiss), and pictures were taken using a cam- 
era (DFC295; Leica) and LAS AF v3.2 software (Leica). For measuring the 
length of dying cells, an aliquot of stained embryos was transferred onto a 
glass slide, and 10 embryogenic structures per line were imaged using 
z-stack and tile-scan functions of a confocal microscope (LSM 780) and a 
20x Plan-Apochromat objective with numerical aperture of 0.8. The length 
of the cells was measured using ZEN 201 1 blue software (Carl Zeiss), and 
>50 cells per structure were analyzed. 

FM4-64, FDA, and DAPI triple staining was performed on 7-d-old 
suspension cultures grown without GF. DAPI and FDA were added to final 
concentration 4 and 2 ug/ ml, respectively. After 15 min, the samples were 
washed twice with — GF medium on a 50-pm nylon mesh. Then, 1 pM FM4-64 
was added, and samples were scanned within 10 min using the sequential 
scanning mode of a confocal (LSM 780) microscope, with excitation of 
405/488/561 nm, emission of 41 0-492/490-587/582-754 nm, and 
a 20x objective with numerical aperture of 0.8. 

Before detection of mRFP-Atg8, samples were treated with 0.5 uM 
ConA (Fluka) for 14 h. If required, samples were additionally incubated in 
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2 pg/ml FDA for 15 min, washed twice with fresh — GF medium, and 
scanned within 15 min after staining. Fluorescence was detected using the 
sequential scanning mode of a confocal microscope (LSM 780), with exci- 
tation of 561 nm and emission of 563-686 nm for RFP and excitation of 
488 nm and emission of 490-587 nm for FDA and an objective of 40x 
water Plan-Apochromat with numerical aperture of 1.2. A 10-pm-thick 
z stack was used to produce a maximal projection image, and mRFP 
puncta were quantified on a randomly selected 10-pm 2 area in ZEN 201 1 
(black edition). The mean intensity of FDA was assessed by selecting at 
least three 10-pm 2 areas per cell using ZEN 2012 blue edition, and au- 
tophagosomes were counted in the selected areas. 

SYTOX orange staining was performed on 7-d-old suspension cul- 
tures grown without GF. SYTOX orange and FDA were added to final con- 
centration 1 pM and 2 pg/ ml, respectively. After 15 min, the samples were 
washed twice with — GF medium on a 50-pm nylon mesh, and samples 
were scanned within 1 0 min using the sequential scanning mode of a con- 
focal microscope (LSM 780), with excitation of 488/561 nm, emission of 
490-587/582-754 nm, and objectives of 10 and 20x Plan-Apochromat 
with numerical apertures of 0.45 and 0.8, respectively. 

TEM 

Samples for TEM were fixed in 2% of paraformaldehyde, 3% glutaralde- 
hyde, and 1 00 mM sodium cacodylate, pH 7.2, for 3 h at room tempera- 
ture followed by overnight fixation at 4°C. Fixed samples were washed 
three times in 100 mM sodium cacodylate, pH 7.2, and postfixed in 1% 
Os0 4 aqueous solution for 1 h at room temperature. Postfixed samples 
were gradually dehydrated in a series of ethanol solutions from 25 to 
100% and infiltrated with LR white resin as previously described in Minina 
et al. (2013a), the resin was cured at 65°C overnight and cut to obtain 
99-nm-thick sections. Sections were examined using a transmission electron 
microscope (CM12; Philips), and images were recorded on film (4489; 
Kodak). Films were then scanned with a resolution of 3,200 pixels/inch 
and processed using Photoshop CS5 extended (Adobe). Raw data were 
rescaled according to magnification used for each image. 

ATP content and oxygen consumption measurement 

7-d-old suspension cultures grown in — GF medium were pelleted on a sterile 
50-pm mesh and washed three times with fresh growth medium to remove 
extracellular ATP. Cells were ground in liquid nitrogen, and ATP was ex- 
tracted by boiling 1 00 mg of the powdered material in 1 50 pi of boiling buf- 
fer (100 mMTris, 4 mM EDTA, and phosphatase inhibitors, pH 7.5). Debris 
from the samples was removed by centrifugation in a table centrifuge for 20 min 
at maximum speed at 4°C. Supernatants were transferred into a 96-well 
plate, and luminescence intensity was measured upon adding the reaction 
buffer (50 mM Tris HCI, 1 mM EDTA, 2 mM MgS0 4 , 7.5 mM DTT, 1 00 pM 
D-luciferin, and 5 pg/ml luciferase from Photinus pyralis) in the plate reader 
(Omega FLUOstar; BMG Labtech). The raw data were then normalized 
using a standard curve as previously described in Minina et al. (201 3a). 

7-d-old suspension cultures incubated in — GF medium were gently 
pelleted by centrifugation and resuspended in 1.6 vol of fresh medium. 
300 pi of the cell suspension was used to measure oxygen consumption 
rate (OCR) in a liquid-phase oxygraph (Oxygraph Plus System; Hansatech 
Instruments) in the absence (OCR normo |) or presence of 10 pM carbonyl cyanide 
m-chlorophenyl hydrazone (OCRd ecoup i e d) or 1 mM KCN (OCR res pi ra tion block). 
Respiratory control ratio was calculated as previously described in Minina 
etal. (2013a). 

Metacaspase activity assay 

Cells from 7-d-old — GF suspension cultures were filtered and ground in liq- 
uid nitrogen. The total protein was extracted using HNCD buffer (50 mM 
Hepes, pH 7.5, 1 50 mM NaCI, 0. 1 % CHAPS, 5 mM DTT, 1 00 pM PMSF, 
and 10 pM pepstatin). Fluorogenic substrate FESR-AMC (final concentra- 
tion of 50 pM; Sundstrom et al., 2009) was added to 7.5 pg of total pro- 
tein in the presence or absence of 50 mM CaCI 2 . Accumulation of 
fluorescence was detected at excitation/emission of 355/460 nm and 
converted into picomoles of AMC released per minute per milligram of 
total protein. Values obtained for reactions without calcium were sub- 
tracted from values obtained in the presence of 50 mM CaCI 2 . Results 
were then normalized relative to control. 

Polyclonal antiserum and immunoblotting 

Antibody production in rabbits was performed according to Canton et al. 
(1996). A total of 100 pg of recombinant mcll-Pa protein mixed 1:1 with 
Freund's complete adjuvant was used for the first immunization. The second 
boost was performed after 40 d, in which incomplete adjuvant was used in- 
stead of the complete version. The antibody was affinity purified using a 
polyvinylidene difluoride membrane containing recombinant metacaspase. 



7-d-old — GF suspension cultures were ground in liquid nitrogen, 
and total protein was extracted in HNCD buffer, mixed 1 :1 with 2x Laemmli 
buffer (Laemmli, 1970), and boiled for 3 min. For Atg8 detection, ex- 
traction was performed in 50 mM Tris, pH 8.0, 150 mM NaCI, 1 mM 
PMSF, and protease cocktail (cOmplete; Roche), and debris from the sam- 
ples was removed by spinning down at 7,000 g for 5 min at 4°C followed 
by filtration through four layers of Miracloth. Cleared supernatants were ul- 
tracentrifuged for 1 h at 100,000 g at 4°C as previously described in 
Chung et al. (2010). Pellets from ultracentrifugation were resuspended in 
extraction buffer containing 0.5% Triton X-100. Proteins were separated 
on a 12% polyacrylamide gel (15% polyacrylamide gel with 6 M urea for 
Atg8) and transferred on a polyvinylidene difluoride membrane. The mem- 
brane was probed with antibodies against NBR1 (kindly provided by 
T. Johansen, University of Tromsa, Tromso, Norway; Svenning et al., 201 1 ), 
Atg8 (Abeam), mRFP (Chromatec), actin (Santa Cruz Biotechnology, Inc.), 
mcll-Pa, or TSN (Sundstrom et al., 2009) diluted 1:1,000 followed by 
anti-rabbit, anti-rat, or anti-mouse HRP conjugates (GE Healthcare) di- 
luted 1 :3,000. The membranes were incubated with ECL Prime kit (GE 
Healthcare), and the signal was imaged using the luminescent image ana- 
lyzer (LAS-3000; Fujifilm). To check protein loading, membranes were 
stained with Coomassie brilliant blue G-250. 

Real-time quantitative PCR (qPCR) 

Total RNA was extracted from 7-d-old suspension cultures grown without 
GF. 1 pg RNA was used per real-time reaction with the Maxima kit (Fer- 
mentas), and 1 /80th of each real-time reaction was used for qPCR analy- 
sis. qPCR reactions were performed using PCR thermal cycler (iQ5; Bio-Rad 
Laboratories) and qPCR kit (DyNAmo Flash SYBR Green; Finnzymes). Prim- 
ers used for qPCR are listed in Table S2. AAC T method was used to quantify 
fold expression of genes of interest in the samples, with two genes, cell divi- 
sion control 2 and phosphoglucomutase, previously reported to show sta- 
ble expression during early embryo development of P. abies (Vestman 
et al., 201 1 ), being used for normalization. 

Statistical analysis 

All statistical analysis was performed in JMP version 10.1 (SAS Institute) 
using either a pairwise Student's t test or comparison with control by 
Dunnett's test. 

Online supplemental material 

Fig. SI shows phylogenetic analysis of Atg5, Atg6, and metacaspase 
homologues from protozoa, fungi, plants, and animals, highlighting the 
position of P. abies homologues. Fig. S2 shows decreased expression of 
ATG5, ATG6, and mcll-Pa in RNAi lines determined by Western blot- 
ting (for mcll-Pa) and real-time qPCR (for all three genes). Fig. S3 supports 
the data presented in Fig. 5 A and shows no effect of the suppression of 
autophagy on mcll-Pa activity using mcll-Pa and TSN (proteolytic target of 
mcll-Pa) Western blot analyses. Table SI presents a list of proteins used 
for phylogenetic analysis shown in Fig. SI . Table S2 shows primers used 
in this study. Online supplemental material is available at http://www.jcb 
.org/cgi/content/full/jcb.201307082/DCl. 
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